Dynamical phase transition in vibrational surface modes 
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We consider the dynamical properties of a simple model of vibrational surface modes. We obtain 
the exact spectrum of surface excitations and discuss their dynamical features. In addition to the 
usually discussed localized and oscillatory regimes we also find a second phase transition where the 
surface mode frequency becomes purely imaginary and describes an overdamped regime. Noticeably, 
this transition has an exact correspondence to the oscillatory - overdamped transition of the standard 
oscillator with a frictional force proportional to velocity. 



I. INTRODUCTION 

In a classical oscillator, dissipation is described consid- 
ering an equation of motion with a friction term propor- 
tional to the velocity [l|. Ignoring microscopic details, 
all environmental effects are summarized phenomenolog- 
ically in the coefficient 770 of this term. We propose here a 
simple and time-reversal invariant model whose analyti- 
cal solution yield dissipation in the thermodynamic limit. 
We consider a variation of the Rubin model [3, yj shown 
in Figfj] a "surface" oscillator (represented as a pendu- 
lum) with natural frequency loq and mass too coupled to 
an ordered and semi-infinite chain of "bulk" oscillators 
whose masses are to. As occurs with a "Brownian bath" 
|^, Ohmic dissipation will require a = m/rriQ <C 1 . 




FIG. 1: Scheme of the model: a simple pendulum (surface 
oscillator) is coupled to the bulk masses. 



The equations of motion are 
d 2 „ . 



dt 2 ' ° 

d2 „ 2 



aw: 



u (t) — au%U\{t) = 0, 



U n (t) - UjI [u n -l(t) + U n +l(t)] = 0, 



(1) 



where to x — y/ ' K/m is the exchange frequency between 
neighbor oscillators and u n (t) denotes the nth oscillator 
displacement from equilibrium. Assuming that the num- 
ber of bulk oscillators is finite, the time-reversal invari- 
ance present in these equations is obvious. We will show 
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below how irreversibility appears in the thermodynamic 
limit where the number bulk oscillators become infinite. 



II. FREQUENCY DOMAIN 

The equations of motion described above can be solved 
by Fourier transforming the displacement: 



u n {t) 



1 

2^ 



and replacing it in Eq.Q: 

(uj 2 1 - R 



due- iu,t u n (u>), 



u(w) = 0. 



(2) 



(3) 



This equation is solved by the eigenfrequencies u>k and 
eigenvectors u(a/fc) whose site components are expressed, 
in bracket notation, in terms of the site versors (n\ as 



u n {wk) = (n\<pk) 
and the dynamical matrix is: 
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We define the Green's function operator: 

ED(w) = (ujH-My 1 , (6) 

whose matrix elements diverge when us coincides with an 
eigenfrequency u>k- 

Consider an initial condition with all masses at equilib- 
rium and an impulsive force mjUj(0)5(t) applied to the 
jth mass, impossing an initial velocity to it. In this case, 
the Green's function provides the time evolution of the 
ith mass displacement: 



«i(«)=-D«(«)«i(0) 



,'\ / due-^D^u)^ ,•(()). (V) 



representing a position-velocity Green's function 
Dij(u>) = D(ui,Uj]{jj), which is related to the position- 
position Green's function: 



D{ui,Uj]uj) = —iojD(ui,Uj;uj). 



(8) 



This gives the displacement of the ith mass when an 
initial displacement Uj(0) is imposed by the force: 



F,{t) = 



%(0) 



[5{t+y)-5{t-\r)\ 



(9) 



at the jth mass, where r -C l/w x - 

In this work, we focus on the surface oscillator through 
Doo(t). If U is the matrix that performs the change of 
basis from sites to diagonal form in M, one has: 

U ik = (i\ (p k ) , (10) 

that is, the fcth eigenmode projection over zth site. In 
this case one can write, 



k 



<Pk) 



(11) 



By performing the analytical continuation uj — ► uj + 
i<5 for Dqo(lj) and taking its imaginary component, one 
obtains the spectral density associated to the surface site. 
In the case where a — > bulk modes and surface modes 
have no projection and oscillations survive undamped. 

Dqq(lu) can be obtained through an infinite order per- 
turbation theory [5| that accounts for the surface mode 
corrections due to the presence of the neighbor oscilla- 
tors. Let us begin with the surface oscillator and a single 
bulk mass. If the two masses are uncoupled (uj x = 0) 
the surface frequency is simply ujq but if u> x ^ 0, Eq.© 
yields 

Ax>(«) = j—. (12) 



Here, ujq is affected by the static correction auj 2 and 
the dynamic correction all (a;) due to presence of the 
other oscillator. 

By taking the limit of the oscillators number to infinite, 
the dynamic correction becomes 



n( w ) = 



2^-1%)' 



(13) 



Because of the presence of an infinite number of oscil- 
lators at the right, a correction on the nth oscillator is 
just the same as that in the (n + l)th oscillator. The so- 
lution of Eq. I|13|) is complex and its imaginary part gives 
the decay rate of a surface excitation with frequency uj. 
Therefore, in the thermodynamic limit, the temporal re- 
currences p (Mesoscopic Echoes jTJj ) disappear. 

III. SURFACE MODE FREQUENCY 

In the case of a surface oscillator with natural fre- 
quency ujq, coupled to an infinite number of bulk oscilla- 
tors, the Green's function is: 



AroM = 



1 



aU(uj) ' 



(14) 



Hence, in the region of continuous modes (|w| < 2uj x ) 
one has 



A)0 M 



1 



1 — § uj 2 — Uq + irj(uj)uj ' 



(15) 



where ujq = u>q/ \fi — a/2 is a first approximation to the 
resonance frequency and 
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2uj r , 



(16) 



describes the dissipation process. Comparing Ea. fl5|) 
with a standard damped oscillator (frictional force pro- 
portional to velocity), it is easy to see that in the broad- 
band limit (a — > 0, oj x — > oo and aoj x = const.) they 
have the same behavior taking t] Q = auj x as friction coef- 
ficient. 

For a dynamic description of the surface oscillator, we 
look for the pole structure of Doo(cu), equivalent to find 
uj such that: 



uj — uj + irj(uj)uj = 0. 



(17) 



If we denote the surface mode frequency as ujr = ^>r. - 
v~in and P as 



P 



<+>Q — U x 

u x \J\ - a 



(18) 



then one has 



ur 



+ Juj% - AujIujI p<-l 



^Juj 2 + 6 [uj 2 _ - i^AujIuj 2 - uj\\ \p\ < 1 



1 + 8 uj 2 _- Juj% - ^uj 2 ujI p > 1 



where uj\ — uj" 2 ± auj 2 and 5 = a/2(l — a). The depen- 
dence of the pole with ujq can be seen in Fig|2 

Here, we show the real and imaginary parts of the pole 
separately and we find three well-defined regimes that 
emphasize the onset of quite different dynamical proper- 
ties. 



IV. DYNAMICAL PHASES 

Localized - extended transition. This transition 
occurs when ujo lays close to the upper band edge. There, 
the real part of the pole makes an excursion into the 
band gap conserving an imaginary component which is 
lost at the cusp zoomed in the inset. This is analogous 
to the virtual states described by Hogreve |8| for elec- 
tronic states which are missed under other approxima- 
tions yj, lllj • They become real normalized states which 
are exponentially localized [llj when the pole touches 
back the band edge. 



a = m/m„ = 0.3 
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FIG. 2: Surface Green's function pole. Above: real part (reso- 
nance). Below: imaginary part (damping). The case lor — loo 
is shown in dotted line. Insets: overdamped - oscillatory tran- 
sition (left), oscillatory - localized (right). 



Considering Eq.JTJ), we can complete our previous anal- 
ysis obtaining the time evolution of surface displacement 
by performing a numerical Fourier transform of the ana- 
lytical expression for D 00 (uj). These results coincide with 
the numerical integration of the equation of motion in the 
Hamilton- Jacobi representation which is performed with 
an ad hoc version of the Trotter-Suzuki method [12j . Set- 
ting the ratio a between masses and Fourier transforming 
for several loq values in the previously mentioned regimes, 
we obtain the displacement pattern shown in FigOU 
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It can be seen in Fig[2]that for large values of surface 
natural frequency (cl>o 3> 2lo x ) we only have real compo- 
nent in lor and the oscillation frequency is approximately 
Wo. In this localized regime, displacements are exponen- 
tially smaller as the oscillators increase their distance to 
the surface. The lack of an imaginary part implies that 
the displacement amplitude, independently of bulk oscil- 
lators, survives indefinitely. In other words, there is no 
energy propagation to bulk oscillators. 

As loo decreases, below a critical frequency lo C 2 = 
lu x (1 + \/l — a), lor becomes complex. Its real part is 
the resonance frequency whereas its imaginary part de- 
scribes the dissipation. This is the extended oscillatory 
regime where surface oscillations have an lor frequency 
and the amplitude of displacement decays with a lifetime 
proportional to 1/"/r. 

Oscillatory - overdamped transition. For loq <C 
lo x the system is in the regime where r)(u>) ~ aw x /(l — 
a/2) and hence Eq. i|17fl yields Ohmic dissipation in which 
the surface kinetic energy decays into bulk modes with- 
out return. This describes a frictional force proportional 
to the velocity. As ujq goes below the critical frequency 
lo c i = lo x (1 — Vl — ck), the pole of the Green's function 
Cjr becomes purely imaginary. This is the overdamped 
regime where no oscillations occur. Notice that 7^ de- 
creases at both sides of the transition, this means that 
a further decrease of loq or increase in the friction 77 also 
implies a decrease in the relaxation rate. All this counter 
intuitive effect arrises from the exact solution of the sim- 
ple mechanical model. 



FIG. 3: Evolution of the surface displacement. White spaces 
denote uo(t) > whereas gray spaces denote Uo(t) < 0. 



Here, the time evolution of the surface displacement 
depends on the initial condition. It is easy to see in this 
picture how gray regions "diverge" near the critical fre- 
quency LV c i so that for smaller coq values uo(t) is always 
positive. On the other hand, the period of the oscilla- 
tion (reflected from the fringe edges) shows a cusp at co C 2 
where it start to increase again. This is consistent with 
the decrease of lor as loq crosses lo C 2 (see inset in Fig|2ll. 



V. CONCLUSIONS 

Dynamical features in dissipation processes have been 
described by a simple classical model that contains the 
essential properties of a surface oscillator interacting with 
bulk vibrational modes. We arrive to results equivalent 
to the phenomenological description of dissipation where 
these effects are summarized in a single term proportional 
to the oscillator velocity. A complete analytical solution 
of the dynamics allowed us to identify a variety of dy- 
namical phases available for the vibrational excitations: 
localized, extended oscillatory and overdamped. Two nu- 
merical methods (FFT and a Trotter-Suzuki algorithm) 
have been employed to obtain the dynamics with equiv- 
alent results. 



The identification of the different dynamical regimes 
becomes very important in nano-physics when cantilevers 
13] involving a relative small number of atoms fail to sat- 
isfy the standard thcrmodynamical approximations. Lo- 
calization, weak diffusion and recurrences become usual 



phenomena in the nanomechanical devices. In a quan- 
tum description of finite size systems, we should consider 
the quantization of oscillation amplitudes. Our formal- 
ism, based on the Green's functions in discrete but open 
systems, seems to be the ideal tool for this purpose. 
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